Inelastic neutron scattering in random binary alloys : an augmented space approach 
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Combining the augmented space representation for phonons with a generalized version of 
Yonezawa-Matsubara diagrammatic technique, we have set up a formalism to seperate the coherent 
and incoherent part of the total intensity of thermal neutron scattering from disordered alloys. This 
is done exacly without taking any recourse to mean- field like approximation (as done previously). 
The formalism includes disorder in masses, force constants and scattering lengths. Implementa- 
tion of the formalism to realistic situations is performed by an augmented space Block recursion 
which calculates entire Green matrix and self energy matrix which in turn is needed to evaluate the 
coherent and incoherent intensities, we apply the formalism to Ni^Pd^ and Ni^Pt^o alloys. Nu- 
merical results on coherent and incoherent scattering cross sections are presented along the highest 
symmetry directions. Finally the incoherent intensities are compared with the CPA and also with 
experiments. 

PACS numbers: 63.50.+X 
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INTRODUCTION 



Phonon excitations in alloys provide a particularly se- 
vere testing arena for any theory of elementary excita- 
tions in disordered systems. The dispersion relations, 
line widths etc. for phonons can be measured directly 
through coherent inelastic neutron scattering. Informa- 
tion about the phonon density of states can be extracted 
from the incoherent inelastic neutron scattering cross sec- 
tions. Over thepast few years, numerous experimental 
studies [l], 0.0a Sill °f the lattice dynamics of disordered 
systems have been carried out providing deep insight into 
the nature of their elementary excitations. However the 
theoritical counter part is still unsatisfactory. The theory 
of scattering of thermal neutrons withperfect crystals has 
been set up on a rigorous basis [HI?], 15 -However, we have 
to make the same statement for disordered alloys with 
more care. This is because the theory of neutron scatter- 
ing in random alloys require two basic inputs : first is the 
formulation of the problem and the second it's actual nu- 
merical implementation in realistic situations. As far as 
the formulation part is concerned, several authors have 
attempted the problem with different approaches. How- 
ever their actual numerical implementation in realistic 
systems still remains sketchy beyond the simple single 
site mean field theories. Quite some time ago, Nowak 
and Dederichs discussed the separation of the coher- 
ent and incoherent parts of the total inelastic scattering 
intensity by using the Yonezawa-Matsubara diagramatic 
technique within the single site coherent potential ap- 
proximation (CPA). According to their approach the in- 
coherent part of the total scattering intensity is the sum 
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of all irreducible diagrams containing only short ranged 
correlations. The coherent part on the other hand may 
be expressed as a product of the usual configuration av- 
eraged Green function and the square of an effective scat- 
tering length, which is itself given by the irreducible di- 
agrams closely related to those for the self energy.lt has 
been known for some time that the single site CPA can- 
not adequately deal with intrinsic off-diagonal disorder of 
the force constants in the problem of phonon excitations 
in random alloys. This was evidenced in the inability 
of the single site CPA to explain experimental lifetime 
data on NiPt Nor can it adequately deal with the 
correlated diagonal and off-diagonal disorder induced by 
the force constant sum rule. Several successful attempts 
have been made to go beyond the CPA. These include, 
among others, approximations based on the augmented 
space formalis m (llOl) : the travelling cluster approxi- 
mation (TCA) lllUld, E3, E3 , the Cluster-CPA [HEi, 
the itinerant cluster approximation (ICPA) an d the 
augmented space recursion (ASR) [la]- 

In this communication we shall tackle a two-fold prob- 
lem : one of formulation and the other of implementation 
in real alloy systems. We shall use a scattering technique 
also based on the augmented space formalism 

013 

to suggest how to separate the coherent and incoherent 
parts of the total inelastic scattering cross-section for a 
disordered binary alloy in a way mirroring the ideas of 
Nowak and Dederichs, but this will be done without tak- 
ing any recourse to mean-field like approximations. For 
implementation in real alloy systems, we shall suggest the 
ASR for the evaluation of scattering cross sections. But 
instead of doing an ordinary recursion, we shall perform 
a Block recursion in order to calculate the off diagonal 
entries of the Green matrix, since the expressions for the 
scattering cross sections in our formalism require the con- 
tribution of off-diagonal Green's function. The approx- 
imation introduced within this formalism will maintain 
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the essential analytic properties of the Green function, 
deal with off-diagonal disorder and the sum rule without 
any further simplifications or assumptions and encom- 
pass environmental effects over an extended neighbour- 
hood. This is the major contribution of this work. 

The rest of the paper is organised as follows. In Sec. 
II, we describe the theory, introducing the augmented 
space representation and it's use in constructing the gen- 
eralized Yonezawa-Matsubara scattering diagrams for the 
averaged green function. In Sec. Ill we derive expressions 
for important physical quantities such as effective scat- 
tering length W e f f (q, w) , inelastic coherent and incoher- 
ent scattering cross sections in terms of the configuration 
averaged Green matrix <C G(q, w) 3> and self energy 
matrix S(q, w) of the system. In Sec. IV, we briefly de- 
scribe the method of Block recursion for the evaluation 
of full Green matrix and Self energy matrix. In Sec. V 
we present our results on Ni^Pd^ and Ni^Pt^ alloys 
as test cases and try to compare them with experimen- 
tal data as far as applicable. Concluding remarks are 
presented in Sec. VI. 



II. FORMALISM 



spanned by the tight-binding basis {|P)|. R,R' specify 
the lattice sites and a, (3 the Cartesian directions, m R is 
the mass of an atom occupying the position R and & R % 
is the force constant tensor. 

We shall be interested in calculating the displacement- 
displacement Green matrix G(R, R',w 2 ) : 

G(R, R',w 2 ) = (R\ (Mw 2 - D) _1 \R') 

Let us now consider a binary alloy A x B y consisting 
of two kinds of atoms A and B of masses ttla and 
tub randomly occupying each lattice sites. We wish 
to calculate the configuration-averaged Green matrix 
<§C G(R, R' ,w 2 ) ^>. We shall use the augmented space 
formalism to do so and indicate the main operational re- 
sults here. For further details we refer the reader to the 
above monograph [2l|]. The first operation is to repre- 
sent the random parts of the secular equation in terms of 
a random set of local variables {n R } which are 1 if the 
site R is occupied by an A atom and if it is occupied 
by B. The probability densities of these variables may be 
written as : 



The augmented space representation for 
phonons 



The augmented space formalism (ASF) has been de- 
scribed in detail in several earlier papers (see [2l||). We 
shall, for the sake of completeness, describe only those 
features which will be necessary for the implementation 
of our ideas in this communication. The basic problem 
in the theory of phonons is to solve a secular equation of 
the form : 

(Mw 2 — D) u(R,w) = 

where u a (R,w) is the fourier transform of u a (R,t), the 
displacement of an atom from its equilibrium position R 
on the lattice, in the direction a at time t. M is the mass 
operator, diagonal in real-space and D is the dynamical 
matrix operator whose tight-binding representations are 

M = ^2 m R <W Pr 

r. 

d = E<« p *+E E n%p RR , 

R R R.'^R. 



along with the sum rule: 



af3 
RR 



E *& 



(i) 



R'^R 



Here Pr is the projection operator \R)(R\ and T RR > is 
the transfer operator \R)(R'\ in the Hilbert space Tt 



Pr(n R ) = x S(n R - 1) + y S(n R ) 

= (-1/tt) Sm(U\ (urI-Nx)- 1 |Tr) (2) 

where x and y are the concentrations of the constituents 
A and B with x + y — 1 . N R is an operator defined on 
the configuration-space 4> R of the variable n R . This is of 
rank 2 and is spanned by the states {|T.r), llfl)}- Wr) is 
a state which indicate that there is an atom A at the site 
R, while || r) is for the atom B. 



Nn 



xp R 



yp l R 



C V T R 



TJ 



2 R 



where p R = \ ] R ){] R \ and p l R = \ |r)(|r | are projection 
operators and t r 1 = \ ] R ){[ R \ and t r ] = \ |r)(|r | are 
transfer operators in the configuration space (j) R . 

In terms of random variables n R , the mass operator 
can be written as : 



r 



5m = mA—m R 
(3) 



According to the augmented space theorem, in order to 
obtain the configuration-average we simply replace the 
random variables n R by the corresponding operators N R 
associated with its probability density, as in Eqn. J5J, 
and take the matrix element of the resulting operator 
between the reference states. For a full mathematical 
proof the reader is referred to ^(| • 
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where 



n R — > N R = x I + (y - x) p L R + Jxy T 1 
Using the above we get, 

M = A(m) 7®/ + B(m) ^V r ®Pr- 

R 

... + F(m) T R ® P R 

R 

= < M > + M' 



U 
R 



A(X) =< X > = (xX A + yX B ) 
B(1Q = (y-x) (X^-Xu) 
F(X) = ^ (Xa - X fl ) 



Similarly the random off-diagonal force constants 
(4) between the sites R and R' can be written as : 



a/3 



<b + (<a + < S -2< b ) n R n R , + (< s - < s ) (n« + n w ) 



(5) 



r 



Let us define the following : 



In augmented space the off-diagonal force constant ma- 
trix becomes an operator : 



a/3 
(1) 
a/3 
(2) 



a/3 



-a/3 
-4B 



'AA "T BB 



2$ 



a/3 
AS 



D 



a/3 
(o//) 



. . . + $g {(y - x) 2 pj, + - *) (pj, T£ +p^ T^) + .xy T£T£} 

E « » ^ ® *W + E ® 



RR' 



The sum rule gives the diagonal element : 



ts?l, = -E{ E 



E E<*>^* 

R I fl'^-R 



D = «D»-£ ^ •■• 

• •• + E 

RR' 

= «D» + D' 



(6) 



The total dynamical matrix in the augmented space is : The boldface operators are 3x3 matrix representations 
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in the three Cartesian directions. 

The augmented space theorem lOj now states 
that the configuration-average of the Green matrix 
<C G(R,R',w 2 ) » may be written as : 

<C G (i?, R', w 2 ) > 

= J G(R,R',w 2 ,{n R }) J] Pr(n R )dn R 

= <{0}®i?| G(w 2 ,{N R })\m®R') 

= ({0}®iZ| (mu^-d)"* |{0}®i?') (7) 

where M and D are the operators which are constructed 
out of M and D by replacing all the random variables 
n R (or n R > ) by the corresponding operators N R (or N R > ) 
as given by Eqn.Q and iJBJ. These are the operators in 
the augmented space = Ti ® The state \R ® {0}) 



is a state in the augmented space, which is the direct 
product of the real-space and the configuration-space 
bases. The configuration-space $ = J}® <p R is of rank 2^ 
for a system of N-lattice sites with binary distribution. 
A basis in this space is denoted by the cardinality se- 
quence {C} — {Ri, i?2, ■ ■ ■ i Rc} which gives us the posi- 
tions where we have a 1 1) configuration. The configura- 
tion {0} refers to a null cardinality sequence i.e. one in 
which we have | "f) at all sites. 

The virtual crystal (VCA) Green matrix is : 

g(R, R', w 2 ) = ({0}®i?| (< M > lo 2 - < D |{0}®i?') 

(8) 

where 

< M >=< m > 7«) I 
Referring back to Equations (Q,© and (0) we get : 



we define 



with 



<G(i?,i?',w 2 ) > = ({0}<»i?| (<M>w 2 -<D>+M'co 2 -D') 1 \{%}®R!) 



= (MSiZlfg^-Dx) |{0}®i?' 



R'jtR 



Pr 



B(m) w 2 + F(m) w 2 T fl n 

n M f„i . n i \ , n (2) ( T n +T n\ , D (3) | i 
+ D (4) f » x T u 4- T u v l )+ D (5) T n T n 



(9) 



(10) 



where 



D (D 

D (2) 
D (3) 
D (4) 
D (5) 



(y - x) ^ 



(i) 



fxy $(i) 



a/3 



(t/-x) 2 $ (2 



= Vxy (y - x) <$>^ 



xy $ 



a/3 
(2) 



B. Generalized Yonezawa-Matsubara scattering 
diagrams for the averaged Green function 

In this section we shall start from Eqn. @ and develop 
a multiple scattering picture based on this. The idea is 
very similar to that of Yonezawa and Matsubara [27J in 
the context of purely diagonal disorder. We shall first 
expand Eqn.® as follows : 
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FIG. 1: The scattering vertices for the averaged Green function 



« G(R, Rf, w 2 ) »= ({0} <g> i?| ( g + g Di g + g D! g Dx g + . . .J |{0} $ it') 

(11) 

I 

Let us discuss very briefly how one generates the scat- The second term yields zero since ({0}<g>i?|Di|{0}<g>i?') = 
tering diagrams. The first term in Equation (jl 1|> gives 0. The third term gives : 

({0} <8> R\ g |{0} <8> i?') = g(i?, i?', w 2 ) (12) 
I 



E E EE <{0}®fl|g|{0}®s') ({01^^'iD! i{c}®^') ... 

S'S" S"'S"" {C} {C} 

. . . ({C} ® S*"| g |{C'} ® 5"') {{C} <g> S'"| f3x |{0} <g> 5"") ({0} ® S"'"| g |{0} ® i?') 



6 



A little algebra yields the following contribution : 

E g(R,S u w 2 ) (Fw 2 ) g(S 1 ,S 2 ,w 2 ) S(S 1 -S 2 ) (Fw 2 ) g(S 2 ,R',w 2 ) 

SiS 2 

E E s(R,Si, w 2 ) Dg } Sa g(S 2 ,S 3 ,w 2 ) SiS, - S 3 ) D^ S4 g(S 4 ,R',w 2 ) 

S1S2 S.3S4 

E E &( R > S *> w2 ) B s!s2 S(S 2 , S 3 , w 2 ) S(S 2 - S 3 ) Dg^ g(S 4 , ^ 2 ) 

E E 8(^.^1. «^) g(S 2 ,S 3 , W 2 ) £(Si - S 4 ) D| 3 S4 g(S 4 ,i?> 2 ) 

S1S2 S3S4 

E E sC^i,™ 2 ) D £s a g(^2,S 3 , W 2 ) <5(S 2 -S 4 ) Dg S4 g(S 4 ,i?> 2 ) 

S1S2 S.3S4 

E E ^, u, 2 ) g(S a , S 3 , w 2 ) <5(S 2 - S 4 )5(Sx - S 3 ) Dg^ g(S 4 , i?', ™ 2 ) 

S1S2 S3S4 

E E 5 i- D £s 2 S(^2, S 3 , u, 2 ) 5(Sx - S 4 )<S(S 2 - S 3 ) Dg S4 g(S 4 , i?', ™ 2 ) 

S1S2 S3S4 



(13) 



■ , 1 



■ l b 



(ii) 



(iv) 



(vi) 



■ , 1 



■ i 



(HI) 



(V) 



(vii) 



FIG. 2: The scattering diagrams for n=2. The top four lines [ 
fig (i) to fig (vii) ] show all possible diagrams. The bottom line 
is a schematic representation of the the topologically distinct 
classes of diagrams 



Referring to Eqns. (|12|l and (|13fl we shall now build up 
the Yonezawa-Matsubara type diagrams. First we shall 
associate scattering vertices with the terms in M' and D'. 
The Fig. ^ shows the seven different type of scattering 
vertices. The dashed lines are associated with the delta 
functions. 

With each factor g we shall associate a propagator repre- 
sented by a horizontal arrow. The connected diagrams to 
order n are then built up by stringing together n+1 prop- 
agators connected by n vertices with all dashed fluctua- 
tion lines connected in pairs. The algebraic terms in Eqn. 
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FIG. 3: The topologically distinct scattering diagrams for 
n=3. 



(|13|l are then represented by the diagram shown in Fig|21 
The Yonezawa-Matsubara diagrams were originally de- 
veloped only for diagonal disorder. The diagrams shown 
in Fig. [SJinvolve off-diagonal scattering terms and the as- 
sociated diagrams are generalized Yonezawa-Matsubara 
diagrams. 

Fig. |2| shows the topologically distinct classes of dia- 
grams for n=3. Note that it involves terms with contri- 
butions from D^ 1 ), D^ 3 ) and D^ 4 ) . These scattering ver- 
tices cannot sit either in the leftmost or in the rightmost 
positions, because one of the associated pseudo-fermion 
Green function line vanishes. 

For n=4, there are various classes of diagrams involving 
D(2) vertices as shown in Fig. 21 In this figure, (i) shows 
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a separable double tent diagram [3l|. The second tent 
goes to renormalize the rightmost phonon Green function 
from g(x,y) to <C G(x,y ^>. (ii) shows a double tent 
non-separable diagram and (iii) a non-separable crossed- 
tent diagram. Of these, the inner tent in the double tent 
diagram (ii) goes on to renormalize the interior Green 
function. As such, the crossed tent diagram (iii) is a 
skeleton diagram of this class. The diagrams (iv)-(vi) are 
a few more examples of skeleton diagrams in this order 
involving other types of vertices. 

If we club together the contribution of all the skele- 
ton diagrams calling this the self-energy, and allow all 
phonon Green functions except the left-most to be renor- 
malized by the separable and non-separable, non-skeleton 
diagrams, we get the Dyson equation : 

<G>= g + gS«G» 

For homogeneous disorder we have shown earlier that we 
have translational symmetry in the full augmented space 
[23| |. We can then take Fourier transform of the above 
equation to get : 



« G(q, E) » = g(q, £)+g(q, E) S(q, E) « G(q, E) » 

The diagrams for the self-energy are skeleton diagrams 
all of which have the structure as shown in Fig. [5] 




FIG. 4: The topologically distinct scattering diagrams for 
n—A. (i) A separable diagram involving only D' 2 ' vertices, 
(ii) A non-separable, non-skeleton diagram involving D^ 2 ' ver- 
tices, (iii) A skeleton diagram involving D^ 2 ' vertices, (iv)- 
(vi) Skeleton diagrams involving various combination of ver- 
tices. 




FIG. 5: Structure of the skeleton diagrams for the self-energy 



III. IMPORTANT PHYSICAL QUANTITIES 

In this section we shall derive results for important 
physical quantities such as the effective scattering length, 
inelastic coherent and incoherent scattering cross sec- 
tions. All of these are derived based on the arguments 
of Nowak and Dederichs for the various kinds of scatter- 
ing diagrams and it turns out that the numerical evalu- 
ation of these quantities require the entire configuration 
averaged Green matrices <C G(q, w) ^ and self energy 
matrices £(q, w) in the reciprocal space representation. 



A. The inelastic neutron scattering cross-section 

The formal expression for the inelastic cross-section for 
the scattering of thermal neutrons from an initial state 
labelled by k to a final state k' with a change of energy, 



E = Hw= J^-{k 2 -k' 2 ) 

and a change of wave- vector q = k — k' + Q , where Q 
is a reciprocal lattice vector is : 
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d 2 a 1 k' 



dVLdE 2NH k 

R W a/3 



EEE «V (w R %mG a R %{w)W R ) n(w) exp {*q • (R - R')} 

(14) 



here : Wr = w R {exp[-(l/2)((q • u R ) 2 ) th ] } : w R is the 
scattering length of the nucleus of the atom sitting at R, 
its equilibrium position, and u^(t) is its deviation from 
equilibrium at the time t. n(w) is the Bose distribu- 
tion function. For a random alloy, Wr, the Debye- Waller 



factor, the atomic mass and the force constants are all 
random variables and dependent on one another via the 
random occupation variables {tir}. Carrying out aver- 
aging over nuclear spins as well as over all the random 
configurations : 



dQ.dE 



1 k' 
2Nh~k 



EEE 

R R' a/3 



Wr G a R %(w) W R , 



)) exp{iq- (R-R')} 



(15) 



Given homogeneity of disorder, we may rewrite the 
above configuration-average as : 



d 2 a 
dVLdE 



1 k' 



2ft k 

av a p 



?V %™ G af3 {^w) n{w) 



Wr G a R p R ,{w) W R ] = G^(R-R',w) 



J2 R g^(R,w) exp{ iq -i?} 

(16) 



Since Wr is a random variable taking two values W a 
and Wb depending on which kind of atom sits at the site 
labeled 'R', So we can write Wr = Wa tir+Wb (1— tir). 
Augmented space theorem then leads to the W factor 
being replaced by an operator in configuration-space as: 



W = A(W)I + B(W)J2 p R®Pr + f ( w )J2 Pr ® T i 

R R. 



-u 

l R 



1 I 
R 



(17) 



A 




V 1 

g g 




w 

1 


w. 



FIG. 6: The scattering vertices related to the fluctuations in 
W 



Since, like the mass, W is a scalar, mode independent 
quantity, A(W) = A(W)J a0 , B(W) = B(W) 5 af} , F(W) 
= F(W) 5 a(} and W = W 5 a/s . 

The scattering vertices arising out of Eqn. (|17fl are shown 
in FigEl 
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The augmented space theorem then gives : 



«W s GBfi-W ff >= (J2® {0}|W (g + gD lg + gD lg Di g 



W|i?'®{0}) 



(18) 



The Fig. [7| shows the scattering diagrams produced 
from Ean. l|18|) for the scattering cross-section. We have 
grouped them into four categories : 

1. In the category (i) are reducible diagrams whose 
end vertices are W or the averaged value «W> 
6 a p- The central decorations are all the diagrams 
we have already seen for the configuration-averaged 
Green function. The contribution of these diagrams 
are therefore : 



4. The last class of diagrams (iv) are irreducible di- 
agrams. Their contribution is also related to the 
self-energy with both the initial and final vertices 
being W 2 rather than T)^ or D 1 ^. We shall dis- 
cuss these diagrams in detail subsequently. 

Because of the disorder renormalization, the W e / f is di- 
agonal neither in real nor mode space and becomes fre- 
quency dependent and complex. The reducible diagrams 
contributes to an expression : 



W < G(w) > W 

2. The second set of diagrams (ii) are also reducible 
diagrams. Inspection of the diagrams immediately 
show us that their contribution may be written as : 

W'(w) < G(id) > W + . . . 
...W <G(w)> w"H 

The contribution of the 'vertex' W'(iy) and W"(u>) 
are closely related to the self-energy, with the ini- 
tial and final vertices being different : W2 rather 
than D2 or D5. We shall discuss this relationship 
subsequently. 

3. The third set of diagrams (iii) are reducible dia- 
grams with contribution : 



(19) 

If we now examine the structure of the self-energy dia- 
grams in Fig. we note that the vertex D2 creates one 
configuration fluctuation at a site if it is an initial vertex 
and annihilates a configuration fluctuation if it is a final 
vertex. On the other hand, the D5 vertex creates two 
configuration fluctuations at two sites if it is an initial 
vertex and annihilates two configuration fluctuations at 
two sites if it is a final vertex. The diagrams for both 
W e / / (w) and the irreducible diagrams in Fig. have 
vertices which create or annihilate only one configura- 
tion fluctuation at both the initial and final vertices. If 
we denote the part of the self-energy contributed by the 
diagrams in the first two lines of Fig. \5\by S(u>), then it 
follows that : 



W'(w) < G(id) > W"(w) 

The contribution of these three sets of reducible 
diagrams may be added together to give : 

W e //H « G(w) » W eff (w) 

where 



W eff (w) =W + W'H 

and 

W o/ /(tfl) = W + W"(» 



W e //(q,uO = Wo + W 2 [A(q)]- 1 S(q, W ) 
Weff(q,w) = Wo + S(q, W )[A(q)]- 1 W 2 (20) 

where 

A(q) = F(m) + D^(q) + D< 5 >(q) 

The expression (|19ll is long-ranged in real-space. Fol- 
lowing the argument of Nowak and Dederichs Q within 
the CPA and Mookerjee and Yussouff 24] in a more gen- 
eral context of a cluster-CPA, we identify contribution of 
the reducible part as the coherent part of the inelastic 
scattering : 
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FIG. 7: The scattering diagrams for the inelastic scattering cross-section 



' d 2 a 
dQ.dE 



coh 



af3 



q a qP 3m 



W e //(q,w) < G(q,tu) > W e //(q,w) n(tu) 



a/3 



(21) 



If we now look back at the irreducible diagrams in class 
(iv) of Fig. |7| we note that the diagrams in the top 
row of (iv) are both short ranged : the leftmost one is to- 
tally diagonal in real-space and the rightmost one has the 
same range as the dynamical matrices, which are reason- 
ably short ranged. The bottom diagram is longer ranged. 
However, this and all long ranged diagrams in this class 
involve crossed-tent diagrams. If we look at the diagrams 
for the self-energy exactly the same kind of structure is 
seen. The only diagrams which can lead to a long-ranged 
self-energy are crossed-tent diagrams like (iii) of Fig. 0] 
These diagrams belong to correlated scattering from con- 
figuration fluctuations at different sites. The contribu- 
tions of such diagrams are dominated by those which are 



short-ranged. Within the single site coherent potential 
approximation (CPA) such diagrams are neglected alto- 
gether and the self-energy (and therefore the irreducible 
diagrams for the cross-section) is diagonal in real-space. 
Beyond the CPA, dominant contributions arise from cor- 
related scattering of the smaller clusters and the contri- 
bution of the irreducible diagrams is also short-ranged 
: being between sites within the smaller clusters. The 
range of the self-energy and the irreducible diagrams are 
therefore as large as the size of the largest cluster whose 
correlated scattering is significant. The Locality princi- 
ple of Heine [25| convinces us that electronic structure 
is insensitive of far-off environment and, although the 
self-energy is not diagonal in real-space (except in the 
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single-site CP approximation), its range is nevertheless 
short. This is behind the reasonable success of the CPA 
in any cases. Again, following the arguments of Nowak 
and Dederichs [3 and Mookerjee and Yussouff |24|. this 
contribution can be related to the incoherent part of the 
inelastic scattering. The incoherent intensity is given by 



d 2 <7 



dVLdE 



where 



^ y r -i a/3 

q a q 13 3to r(q, w) n(w) 



2h k 



a/3 



(22) 



finite scatterers attached to perfect leads. We shall bor- 
row their ideas and set up a block recursion in the space 
of vibrational modes (rather than the lead space, as in 
Godin and Haydock's work) in order to obtain the Green 
matrices directly. 

The recursion method essentially starts from a denu- 
merably infinite basis and changes the basis to one in 
which the dynamical matrix (or the Hamiltonian, in elec- 
tronic problems) is tri-diagonal. In the block recursion 
we start from a matrix basis of the form : where 
J is the discrete labelling of the augmented space states 
and the a, (3 labels Cartesian directions (i.e. the modes 
of vibration) . The inner product of such basis is defined 
by: 



r(q,«,) = W 2 [A(q)]- 1 H(q,w) [A(q)]- 1 W 2 (23) 



IV. BLOCK RECURSION AND THE 
SELF-ENERGY MATRIX 



(*(»>,*<«>) = *wt 



(™)t fls(m) 
,7 V J,P'P 



i\jnm 

— iV a/3 



J p' 



For a real-space calculations on a lattice with Z- nearest 
neighbours, we start the recursion with : 



The augmented space recursion (ASR) has been pro- 
posed earlier by us [Til |26| as technique for the incor- 
poration of the effects of configuration fluctuations for 
random substitutionally disordered alloys. This can be 
achieved without the usual problems of violation of the 
herglotz analytic properties [Tg of the approximated con- 
figuration averaged Green functions for phonon problems 
in realistic random alloys. Although our initial focus 
was on spectral functions and complex dispersion rela- 
tions and lifetimes, in this communication we propose to 
study inealstic neutron scattering cross-sections. These 
calculations require the full Green matrices and not only 
their diagonal elements. We propose here the use of a 
generalization of the recursion method of Haydock et.al. 
|28( . The block recursion technique had been introduced 
earlier by Godin and Haydock |2JJ, |3ij in the very dif- 
ferent context for obtaining the scattering S-matrix for 



(i) 

J,ctp 



Uj$5j tl + U$5j t2!+1 



(2) 



while for a reciprocal space calculation we start with : 



a i _ TT W : ,77 2 a 

$ J,a/3 - U a (3 °JA + U a p °J,2 



where 
U 



V x, _ ag ■ , r . rA2) _ F( m -^ 2 ) c 



The remaining terms in the basis are recursively ob- 
tained from : 



Z^^J,a/3'^f3'l3 ~ 2s 2s Uja < J & 2^* J,a/3> A P' 13 

13' ,/' /3' (3' 

®J,a0> B p>f3 ~ 2s2s H J^J'P'®J',[3't3- 2^®J, a f}' A p'f}- 2-r*W B P'P 



P' 



J' p> 



(24) 



where, H = M 1 / 2 D M 1 / 2 with M and D given by Orthogonalization of the basis gives : 
Eqns.lQJ and @J. 



EEEE = E ^ A P 



j p' j' p" 



p' 



12 



In matrix notation, where matrices are in the vibrational 
mode (a/3) space : 

A (n) = fanny 1 J^J^ ^ H./J' ( 25 ) 
J J' 

Next, we note that we had started with a orthogonal 
basis set of rank J max x a max . The above procedure 
merely gives Jmax basis sets. We still have orthogonality 

(n) 

conditions among the various columns of ^j^g- in order 



to impose these conditions, consider 

= E E nja,j>fp*%f, - E ^ 41 ■■■ 

J' (3> P> 

Construct three column vectors tpf^ out of the three 
columns of ^!j ta p and set about to Gram-Schimidt or- 
thonormalizing the set : 



a J a J 



^ = E^» =* ^ = E^V^ (*<»») ; s^E^V^-E^ 

(26) 



km 

k—1 aj exj k—1 



where m stands for the total number of vibrational 
modes. 

We may now construct ^ j^g out of (f>j a and note that 

Bkm is indeed the matrix B^ n+1 ^ we are looking for. 

The Eans. l24(l - H26|l show that we may calculate the 
matrices {A^ n \ B( n+1 )} recursively, noting that B^ 1 ) = I 
and B^ ) = 0. In this new basis, the Hamiltonian is 
block tri-diagonal and the Green matrix can be written 
as follows : 



B t Q(n+1) B 



for n > N 



A judicious choice of S (0.001) and N (5000) gives a 
smooth density of states from the diagonal part of the 
Green matrix. The self-energy follows from the Dyson 
equation : 

S = g- 1 - G 1 



Q(n) 



w 



2 I A (n) - B ( " +1)t G (n+1) B (n+1) 



«G» = G (1) 



(27) 



The terminator which replaces the asymptotic part of the 
matrix continued fraction is that which is used by Godin 
and Haydock |3(|. We calculate the matrix coefficients 
upto a n = Nq and approximate at coefficients > Nq by 
A and B. We then write for a N > N Q : 



G (N) = 



(w 2 - iS)I 



and then iterate 



RESULTS AND DISCUSSION 



In the next two subsection, we explore the relative im- 
portance of mass and force constant disorder in the in- 
elastic neutron scattering for two specific random alloys 
Ni^Pd^ and Ni^oPt^o- In a recent publication 0], we 
have already studied the advantages of Augmented space 
recursion (ASR) over the simple CPA for understanding 
the dispersion and life time of phonons in random binary 
alloys. The present work is an extension of that work 
from the implementation point of view, since now we 
need to apply the Block recursion technique to calculate 
the full Green matrix. 
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Frequency v(THz) 



FIG. 8: The coherent scattering cross section in different di- 
rections for NissPd^s. In each of the different directions, the 
various curves indicate the cross sections for various £ values 
starting from the lowest value to the edge of the Brillouin 
zone. The y-axis is in an arbitrary scale with heights scaled 
to the maximum height. Different curves for different £ values 
are shifted along the x-axis in order to facilitate vision. 



Ni55Pd45 alloy : Strong mass and weak force 
constant disorder 



The mass disorder in this alloy is much larger than 
the force constant disorder. We refer the reader to a 
previous article 18] by us, for some of the basic proper- 
ties of fee Ni and Pd, which is relevant for our present 
calculation. This particular alloy has already been stud- 
ied experimentally by Kamitakahara and Brockhouse Q 
using Inelastic neutron scattering. The properties as- 
sociated with the neutron scattering of Ni and Pd are 



reasonably favorable. The incoherent scattering cross 
section of Ni is fairly high, which was at first thought 
to be a potentially serious disadvantage, but in prac- 
tice, this did not turn out to be much of a problem, it 
has been found that the scattered neutron distributions 
are always dominated by the coherent scattering even for 
high frequency made with large widths. Experimental in- 
vestigation shows that significant difference between the 
coherent scattering lengths of Ni and Pd (The coherent 
scattering length for Ni is 1.03 while that of Pd is 0.6 in 
units of 10~ 12 cm ) produces additional incoherent scat- 
tering in Ni^Pd^, but this is much smaller than the 
incoherence produced by the Ni-atom itself in the alloy. 

In Fig.JSJl, we display the inelastic coherent scatter- 
ing cross sections [calculated from Eqn. I|21(l ] obtained 
from our Block recursion calculation along the highest 
symmetry directions ([COO], [CCO], [CCC] , C = |<z1/|<Zmaa|)- 
For a particular direction, different curves indicate the 
cross sections at various ^-points starting from the low- 
est value (C = 0) to the edge of the Brillouin zone (( = 1 
in units of 27r/a). The first thing to note is that the scat- 
tering cross sections are often asymmetric near the res- 
onances. This property was also reflected in the phonon 
line shapes shown in the previous communication [Lsl ] by 
us. However the amount of asymmetries in the cross sec- 
tion is more than that in the usual lorentzian phonon 
line. That should be obvious because if we see the pa- 
per of Nowak and Dederichs Q ; there they have derived 
an expression for the coherent scattering cross section in 
the single site CPA framework, this expression contains 
in addition to the usual lorentzian phonon line contri- 
bution [ obtained from 3m <C G(q,w) ^> ], a second 
term which will have zero contribution only if the scat- 
tering length do not fluctuate ( which is not the case in 
our formulation ) . They have also argued that this extra 
term leads to an asymmetry (rather small contribution) 
of the phonon line. One can also notice that the nature 
of asymmetry in the cross section is not the same as in 
the line shapes. This is due to the contribution of off- 
diagonal elements of the Green's matrix and self energy 
matrix in the coherent scattering cross section [ as obvi- 
ous from expression l|21l) ] . This difference in the nature 
of asymmetry is more pronounced in the [CCO] and [CCC] 
symmetry directions, because the Green's matrix [ and 
self energy matrix ] comes out to be completely diagonal 
in the [COO] direction. The occurence of such a structure 
of the cross section may also be due to the calculation 
in the mixed mode frame work. One should notice from 
the Block recursion technique described in section (IV) 
that, unlike the ordinary recursion where one extracts 
results for specific modes, the Block recursion requires 
a mixed mode starting state and hence does the calcu- 
lation in that frame work to evaluate the entire Green 
matrix. The asymmetries can be described as a tendency 
of more scattering to occur near the resonance frequen- 
cies. It is important to note that the coherent scatteirng 
cross sections have a pronounced q-dependence in all the 
three symmetry directions. Because of the short range 
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Frequency v(THz) 

FIG. 9: The incoherent scattering cross sections in different 
directions for Ni^Pd^ alloy with <&Ni_pd = 0.7 <&Ni-Ni ■ 
In each of the different directions, the various curves indicate 
the cross sections for various £ values starting from the lowest 
value to the edge of the Brillouin zone. 



properties, the self energy ( and [A(q)] _1 ) depends only 
rather weakly on q and does not show any strong struc- 
ture as a function of the same. The same applies for the 
effective scattering length W e //. Thus the only strong 
q-dependence in the coherent cross section arises from 
the average Green's matrix <C G(q, w) ^S> which is a 
long range matrix due to it's dependence on reducible 
diagrams. 

In Fig.Q we display the incoherent scattering cross 
sections [ calculated from Ean. (|22[l ] along the highest 
symmetry directions. In a particular direction, the dif- 
ferent curves indicate the cross sections for various £- 
points starting from £ = to (, = 1 (in units of lixja). 
A look at fig.® immediately shows that the incoherent 
cross sections are very weakly dependent on q. It is the 




2 4 6 8 10 

Frequency v(THz) 

FIG. 10: The incoherent scattering cross sections in different 
directions for Ni^Pd^ alloy with &Ni-Pd = 0.9 <&Ni-Ni • 
In each of the different directions, the various curves indicate 
the cross sections for various £ values starting from the lowest 
value to the edge of the Brillouin zone. 



f ac tor in Ean. (l22|) which weights up the cross 
sections as we go on increasing q-points. These results 
are in accordance with the arguments of Nowak and Ded- 
erichs @ and Yussouf and Mookerjee |24j. They have also 
mentioned in their paper that, the weak q-dependence of 
incoherent scattering cross section arises because of it's 
strong similarities with the self energy diagram which is 
itself a short range matrix due to it's dependence on ir- 
reducible diagrams and hence vary rather weakly with 
q. Kamitakahara and Brockhouse also found a similar 
qualitative features for the coherent and incoherent cross 
sections in their Inelastic neutron scattering measurment. 

Fig. (|1U|> shows the same Incoherent scattering cross 
sections for Ni^Pd^ alloy but with a different 
parametrization of force constants. Here we keep 
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a/3 

Ni-Ni 



and Qp^-pd the same as those of the pure 



materials (same as in Fig|§J) and reduced the & c Ni _ Pd 



s q/3 



below the &Ni-m 



^Ni-Pd 



= 0.9 $ 



a/3 

Ni-Pd 



7 <J> Q/3 
u -' * Ni-NV 



by an a/3-independent factor as 
. However in Fig.© it was 
One can easily see the varia- 



a.0 

Ni-Ni 



tion in shape of the cross section as the input parameter 
varies. It is because of this reason a prior information 
about the species dependence of the force constants is 
very important. Our approach here made no attempt to 
obtain the input parameters themselves from first prin- 
ciples, but rather resorted, as others did earlier, because 
the aim of this work is to establish the augmented space 
block recursion as a computationally fast and accurate 
method for the cross section calculation in context of 
phonon excitations in random alloys. Our future en- 
deavor would be to rectify this, and attempt to obtain 
the dynamical matrix itself from more microscopic the- 
ories, so that we can have a unique shape of the cross 
sections unlike the present case. 



NisoPtso alloy 



Strong mass and force constant 
disorder 



In this section we shall apply our formulation to NiPt 
alloys where both kinds of disorders are predominant. 
The mass ratio mpt/m^i is 3.3 (quite large compared 
to that in NiPd system) and the force constants of Pt 
are on an average 55% larger than those in Ni. For a 
list of general properties of fee Ni and Pt, we refer the 
reader to article jig. Tsunoda et.al. p| investigated 
Ni-\_- x Pt x by Inelastic neutron scattering and compared 
their observations with the CPA. Here for illustration, 
we have considered x = 0.5 only because that makes it a 
concentrated alloy and the failure of the CPA was, qual- 
itatively more prominent at this concentration. Because 
of the large mass and force constant differences, the ef- 
fect of disorder in NiPt alloy is dramatic, such as the 
appearence of sharp discontinuities (or split bands) ob- 
served in the dispersion and line width |l8j . The present 
theoritical investigation also found this kind of peculiar 
behaviour in the inelastic scattering cross sections. 

In Fig. (|llf> . we display the inelastic coherent scattering 
cross sections along the highest symmetry direction. As 
before, In a particular direction the different curves indi- 
cate the cross sections for various ^-values. For the sake 
of simplicity, we have used the same parametrization of 
masses and force constants as used in our previous paper 

m 

In Ni^oPtm ahoy, the coherent scattering cross sec- 
tions show few extra features. Even in the [COO] direction, 
the cross section becomes well seperated double peaks 
along with weakly defined peak in between in the region 
from £ = 0.68 to the zone boundary. The occurence of 
such a weakly defined peak is due to the inclusion of force 




Frequency v(THz) 

FIG. 11: The coherent scattering cross section in different 
directions for NisoPtso- In each of the different directions, the 
various curves indicate the cross sections for various £ values 
starting from the lowest value to the edge of the Brillouin 
zone. 



constant disorder explicitly in our formulation. It is also 
clear from Fig. 1)11(1 that there exists no appreciable peak 
intensity below 3.7 THz for £ > 0.68 in all the three sym- 
metry directions. Tsunoda also found the same structure 
below 3.5THz for C > 0.7. However CPA predicted the 
lower frequency peak to exist for all the C-values. 

For smaller £- values, the lower frequency peaks are 
sharper than the high frequency ones, but the intensity of 
former decreases significantly with increasing wave vector 
however the latter gets sharper, in all the three symmetry 
directions. The phonon peaks are well defined for smaller 
(C < 0.38) and higher ^-values, but no well defined peaks 
were observed for intermediate values of C, presumably 
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Frequency v(THz) 

FIG. 12: The incoherent scattering cross section in different 
directions for Ni$oPt$o. In each of the different directions, the 
various curves indicate the cross sections for various £ values 
starting from the lowest value to the edge of the Brillouin 
zone. 



due to extreme line broadening. This kind of qualitative 
features has also been observed experimentally by Tsun- 
oda et. al. for longitudnal branches. This feature is more 
transparent in the [CfO] direction [ fig. QUE ] where the low 
frequency peak gets broadened and becomes more asym- 
metric in the region between £ = 0.4 to £ = 0.78. 

The incoherent scattering cross sections (given by 
Ean l22[l for Ni^Pt^o alloy along the highest symmetry 
directions are shown in Fig ll2l The weak q-dependence 
of the cross section is obvious from the figure. The in- 
tensities for various q's in a particular direction have 
nearly an approximate q 2 dependence because of the fac- 
tor 'q a qP' in Ean.p5|l. 

In Fig. i|13fl . we compare our results for the incoher- 
ent scattering cross section with those of the CPA and 



the experiment 1]. Here the left, middle and right panel 
displays the augmented space block recursion (ASBR) 
result, experimental curve and the CPA result respec- 
tively. In the CPA result we can observe a dip at the fre- 
quency corresponding to the phonon band gap observed 
in the dispersion curves. This suggests a split band be- 
haviour which clearly seperates the Pt-contribution in 
the low frequency region from the Ni-contribution in the 
high frequency region, because the low frequency region 
is dominated by the Pt-atom (heavier atom) having much 
lower incoherent scattering length than Ni [ the incoher- 
ent scattering length for Pt is 0.1 while that of Ni is 4.5 
]. This kind of spurious gap however is not observed in 
the recursion result, because the CPA results are based 
on the mass fluctuations alone, ignoring the off-diagonal 
and environmental disorder arising out of the dynamical 
matrix. On the other hand, by incorporating the force 
constant disorder as is done in the bolck recursion, we get 
rid of this spurious gap and obtain rather a good agree- 
ment with the experimental results. The overall qualita- 
tive behaviour is similar. In addition the phonon band 
edges in the recursion results are very close to the ex- 
perimental ones. The recursion finds the right band-edge 
at ~ 7.91THz, Tsunoda finds this band-edge experimen- 
tally at ~ 7.93THz while the CPA gives a rather higher- 
value of ~ 8.267 THz. 



VI. CONCLUSIONS 

We have presented a straightforward and tractable for- 
mulation for the seperation of total intensity of thermal 
neutron scattering from disordered alloys into a coherent 
and an incoherent part. The use of the augmented space 
to keep track of the configuration of the system has made 
the formalism simple yet powerful. In essence, the split- 
ting is identical to that introduced by Nowak and Ded- 
crichs 9] within a Yonezawa-Matsubara diagram tech- 
nique except that it has been done exactly without taking 
any recourse to mean-field like approximation. Unlike the 
method proposed by Yussouff an Mookerjee (24|, where 
the diagram technique was exceedingly difficult to gener- 
alize into even a small cluster CPA, the augmented space 
block recursion proved to be simpler to apply the for- 
malism on realistic random alloys. The technique takes 
into account fluctuations in masses, force constants and 
scattering lengths of the individual nuclei. The environ- 
mental disorder arising out of the force constant sum rule 
has also been incorporated before averaging. The ap- 
proximation involving termination of matrix continued 
fraction expansion of the Green matrix retains the essen- 
tial Herglotz analytic properties of the diagonal Green's 
function. We have applied the method to NiPd and NiPt 
alloys. In Ni^Pd^ 1 we have demonstrated that mass 
disorder plays the prominent role. In addition our coher- 
ent scattering cross section enable us to understand the 
effect of small contribution of the off-diagonal elements 
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FIG. 13: The incoherent scattering cross section for VisoPtso. The left panel, middle panel and right panel display the block 
recursion result, experimental curve and the CPA result respectively. 



of Green matrix. The results on Ni^oPtso alloy however 
demonstrate the prominence of force disorder even in a 
case where the mass ratio is ~ 3. Our results agree well 
both with the coherent and the incoherent scattering ex- 
periments, where as the CPA fails both qualitatively and 
quantitatively. We propose the technique as a computa- 
tionally fast and efficient method for the study of inelastic 
neutron scattering in disordered systems. Our approach 



here had no prior information about the species depen- 
dence of the force constants, but rather choose a set of 
force constants intutively as others did earlier. A bet- 
ter understanding of the role of disorder in the lattice 
dynamics of random alloys could be achieved with prior 
information about the force constants. These could be 
obtained from more microscopic theories ( e.g. the first 
principles calculation on a set of ordered alloys.) 
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